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We describe results of measurements of the orientational motion of glass microrods in a microchannel 
flow, following the orientational motion of particles with different shapes. We determine how the 
orientational dynamics depends on the shape of the particle and on its initial orientation. We find 
that the dynamics depends so sensitively on the degree to which axisymmetry is broken that it is 
difficult to find particles that are sufficiently axisymmetric so that they exhibit periodic tumbling 
(‘Jeffery orbits’). The results of our measurements confirm earlier theoretical analysis predicting 
sensitive dependence on particle shape and its initial orientation. Our results illustrate the different 
types of orientational dynamics for asymmetric particles predicted by theory. 


I. INTRODUCTION 

We study experimentally the orientational dynamics of neutrally buoyant non-axisymmetric particles 
suspended in a viscous shear flow. The rotation of axisymmetric particles in a shear flow has been studied 
in several experiments (we give a brief account below). The case of non-axisymmetric particles, by contrast, 
has received little attention experimentally. This is surprising because it is known in theory^’^ that the 
orientational dynamics can be very sensitive to small deviations from the axisymmetric limit. How sensitively 
the dynamics is affected by slight breaking of axisymmetry depends npon the orientational trajectory which 
in turn is determined by the initial orientation. 

In order to verify this theoretical prediction experimentally it is necessary to use particles with well- 
defined shapes, to make sure that inertial effects and thermal noise are negligible, and to compare different 
orientational trajectories of the same particle. 

In this paper we describe experimental observations of the orientational motion of micron-sized glass par¬ 
ticles suspended in a pressure-driven micro-channel flow. The particles have different shapes: axisymmetric, 
slightly non-axisymmetric, and substantially non-axisymmetric (the latter are strongly triaxial particles 
made for the purpose of this experiment by joining two micron-sized glass rods). We verify that inertial and 
thermal torques have negligble effects by showing that the orientational dynamics is invariant under reversal 
of the pressure. An optical trap allows us to manipulate the same particle into different initial orientations, 
to study different orientational trajectories of the same particle. The results of our measurements confirm 
the predictions of Refs. 1,2, show the sensitive dependence of the orientational dynamics upon particle shape 
and initial orientation, and qualitatively illustrate the different types of orientational dynamics computed in 
Refs. 1,2. We have not attempted to quantititatively compare individual experimental particle trajectories 
with theory, to compare individual trajectories requires to compute the resistance tensors for the actual 
shape of the particles used in the experiments (the theory of Refs. 1,2 is formulated for ellipsoidal particles). 
This is beyond the scope of the present paper. Besides the resulting trajectories are of little general interest, 
they depend very sensitively on the particular shape of the particle in question. 

The remainder of this Introduction briefly introduces the relevant theory as well as the wider context of 
this work. The orientational dynamics of a particle suspended in a viscous flow is determined by resistance 
tensors that relate the local flow velocity and its gradients to the torque acting on the particle^’^. A given 
particle shape corresponds to a set of resistance tensors. Their elements are computed by solving Stokes’ 
equations in appropriate geometries^’^. The case of an ellipsoid in a viscous shear flow was first solved 
by Jeffery^. In this limit the equation of motion is given by the condition that the hydrodynamic torque 
vanishes at every instant. For particles that possess three orthogonal mirror planes®, particle shape enters 
the orientational equation of motion through two parameters, A and K. For an ellipsoid with half-axes a, b 
and c, for example, A = (A^ —1)/(A^-|-1) and K = — 1)/-\-l) with aspect ratios X = a/c and k = hjc. 


1 These authors contributed equally to this work. 
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The case K = 0 corresponds to an axisymmetric particle. For axisymmetric particles the orientational 
dynamics is exactly solvable^. When |A| < 1 (so that 0 < A < oo) there are infinitely many degenerate 
periodic orbits, the so-called ‘Jeffery orbits’. This is a consequence of the fact that the dynamical system 
has a conserved quantity: the ‘orbit constant’. The value of the orbit constant is determined by the 
initial orientation of the particle. The orientational motion of an axisymmetric particle in a simple shear 
is sometimes referred to as ‘tumbling’, the particle spends a long time aligned with the flow direction, and 
it periodically changes orientation by 180 degrees. Different Jeffery orbits differ in the functional form of 
these periodic ‘flips’. The degeneracy for iF = 0 is particular to the simple shear flow, and it means that 
small perturbations can have a large effect. 

Inertial forces, for example, are neglected in Jeffery’s theory. These forces induce ‘orbit drift’ into a 
final stable orbit. This was already suggested by Jeffery^, and was discussed in many experimental papers 
starting with Taylor^. See also Ref. 8. The corresponding theory is discussed by Subramanian and Koch® 
and by Einarsson et 

Small particles may be affected by thermal noise so that the orbit constant performs a random walk giving 
rise to a statistical distribution of orientations. This mechanism forms the theoretical basis for understanding 
the rheology of dilute suspensions^®’^^. 

A third possibility, the topic of this work, is that the particle is not perfectly axisymmetric. This leads to 
a more complicated orientational equation of motion, also derived by Jeffery^. Some numerical examples of 
its solutions were reported by Gierszewski and Chaffey who found that the motion of a non-axisymmetric 
particle in a simple shear flow is qualitatively different from that of an axisymmetric particle. Hinch and 
Leal ^ analysed the structure of the solutions to the equation of motion. They found that for short times 
a nearly axisymmetric ellipsoidal particle approximately follows a Jeffery orbit, but on longer time scales 
the ‘orbit constant’ does not remain constant. It oscillates, giving rise to ‘doubly periodic’ tumbling: time 
series of the components of the unit vector aligned with the major axis of the particle show two distinct 
periods. Subsequently Yarin et al.'^ inferred from numerical experiments and analytical calculations that 
ellipsoidal particles may tumble periodically, quasi-periodically, or in a chaotic fashion - depending on the 
particle shape and on its initial orientation. The term ‘quasi-periodic’ refers to doubly periodic motion 
with incommensurable periods. Our experimental results support these theoretical predictions: our analysis 
demonstrates that the tumbling may indeed be periodic, doubly periodic, or possibly chaotic, depending 
on particle shape and initial orientation. Our results are in good qualitative agreement with theoretical 
predictions. 

This work considers the orientational motion of small neutrally buoyant particles in a time-independent 
viscous shear flow. This is a special but important case. It is of theoretical interest because of its degeneracy 
and sensitivity to small perturbations, and it is of practical interest because it fundamentally relates to 
theories and experiments concerning the rheology of dilute suspensions. Theories are commonly formulated 
in terms of Jeffery’s equation^®’^®. Recently there has been a surge of interest in describing the tumbling of 
non-spherical particles in turbulentand other complex flows®^^®^. Since it is difficult to solve the coupled 
particle-flow problem most theoretical and numerical studies rely on Jeffery’s equation as an approximation 
to the orientational dynamics. Some exceptions are described in Refs. 26-28. 

This article is organised as follows. In Section II we enumerate previous experimental efforts to validate 
Jeffery’s equations. In Section III we describe the experimental setup. Section IV contains our experimental 
results. These results are discussed in Section V, and we conclude in Section VI. 


II. PREVIOUS EXPERIMENTAL WORK 

In this section we give a brief account of earlier experiments observing the orientational dynamics of single 
particles in shear flows. 

Taylor ^ immersed millimeter-sized aluminum spheroids in sodium silicate between two concentric rotating 
cylinders approximately 10 mm apart. In his brief report he asserts that the tumbling of the spheroids is 
in qualitative agreement with Jeffery’s predictions, but that the orientational dynamics drifts after many 
(order of 100) particle rotations. 

In a related study Eirich et observed the orientations of glass rods and silk fibres in a Taylor-Couette 
device. The particles had diameters between 10 and 50 pm and aspect ratios between 5 and 100. No 
quantitative data on the orientational dynamics was measured, but they observed that the particles tend to 
align with the flow direction or with the vorticity direction. 

Binder^® studied fibers of many different aspect ratios suspended in glycerine. He employed a similar 
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device with two concentric cylinders and found, as Taylor, that the orientational dynamics slowly drifts. 

Mason and co-workers have studied the orientational dynamics of small particles in shear flows during 
two decades^^^^®. Initially Trevelyan and Mason used a setup of two concentric cylinders rotating in 
opposite directions, making it possible to study a single particle over an extended period of time. The 
gap between the cylinders was approximately 10 mm, the suspending liquid was white corn syrup, and the 
particles were 9.5 p.m diameter glass fibers cut to different lengths. By observing the particle orientation 
in a plane orthogonal to vorticity, Trevelyan and Mason found fairly good quantitative agreement of 
their experimental results with Jeffery’s theory for one particle rotation (Fig. 7 in Ref. 31). However, 
for longer time series (up to 30 revolutions) their results were inconclusive: sometimes orbit drift was 
observed, sometimes not, and sometimes the change of orbit appeared seemingly erratic. In order to compare 
quantitatively with Jeffery’s equations, valid for spheroidal particles, Trevelyan and Mason fitted the value 
of A to measurements of the tumbling period, yielding in their words an ‘effective aspect ratio’ for cylinders. 
Bretherton® later showed that this procedure is consistent. Mason and Manley®^ extended the experiment to 
hundreds of particle rotations (Fig. 2 in Ref. 32), but the observed orbit drift was apparently erratic. Mason 
and Manley mention convective currents as a possible cause for the observed drift, but no conclusions could 
be drawn concerning the single-particle dynamics. Bartok and Mason®® used a similar device consisting 
of concentric cylinders, with a camera-equipped microscope observing along the vorticity axis, allowing to 
very precisely measure the tumbling behavior of high aspect ratio (A « 45) acrylic (‘Orion’) fibres. The 
experimental results were found to agree quantitatively with Jeffery’s predictions for one particle rotation 
(Fig. 5 in Ref. 33). However, no data on the orbit drift, if any, was presented. Goldsmith and Mason®^ 
performed the first quantitative measurements of the rotations of disks. They used the same coaxial-cylinder 
setup described above, with silicone oil for the suspending liquid. The disks were fabricated by heating and 
compressing polystyrene spheres. The disk diameters were 400-850 |um, and their aspect ratios ranged 
from A = 1/20 to A = 1/4. Goldsmith and Mason showed that Jeffery’s theory quantitatively predicts 
the orientational motion of axisymmetric disks, and that the orbit remained constant over 120 particle 
rotations (Table H in Ref. 34). In a sequel Goldsmith and Mason®® described measurements of the motion 
of particles suspended in a liquid flowing through a circular tube. The tube diameters were 2-8 mm, the 
flow was pressure-driven by a syringe pump, and observations were recorded through a microscope traveling 
along the tube. Measurements were made on an assortment of particles of sizes wO.l mm with aspect ratios 
ranging from A = 1/20 (disks) to A = 100 (fibres). Goldsmith and Mason concluded that the dynamics along 
Jeffery orbits compares well with theory for short times (Fig. 5 in Ref. 35). However, they found measurable 
orbit drift after a single particle rotation, for both a rod and a disk, see Fig. 7 in the same paper. The drift 
was attributed to particle asymmetry. Anczurowski and Mason®® fabricated prolate spheroidal particles by 
polymerising an electrostatically deformed droplet, and showed that Jeffery’s theory holds quantitatively for 
one particle rotation given the true aspect ratio A. They used the same concentric-cylinder device described 
above. 

Harris et performed experimental measurements on non-axisymmetric particles {K ^ 0). They 
used an apparatus with counter-rotating cylinders with a gap of 27 mm, which was filled with a glucose 
solution. The particles were machined from a composite material (‘Tufnol’) into cuboids of 1.75 mmxl.28 mm 
cross-section and 2.5-9.5 mm in length. They measured the unknown elements of the resistance tensors by 
observing simple rotations around each of the principal axes. With these numerical values of the resistance 
tensors they compared orientational trajectories of the cuboids to numerical solutions of Jeffery’s equations 
and found reasonable quantitative agreement for two particle revolutions (Fig. 9 in their paper). 

Stover and Gohen ®® investigated the effect of a wall on fibre motion using a pressure-driven flow of corn 
syrup through a rectangular channel. The fibres had cross-sectional diameters of 50 qm and lengths 600 qm, 
resulting in an aspect ratio of 12. They found that the orientational dynamics are in quantitative agreement 
with Jeffery’s theory for two particle rotations when the particle is at least one particle length away from 
the wall. 

Kaya and Koser ®® observed E. coli bacteria advected in a microfluidic channel of rectangular cross section. 
They found that the orientational motion of the bacteria approximately follows Jeffery orbits. 

Einarsson et al. described examples of orientational trajectories of polymer microrods in a microchannel. 
They found that the trajectories of some particles were periodic, admitting comparison to Jeffery’s theory. 
Other trajectories were seen to be aperiodic (Fig. 8 of Ref. 27). The authors suggested that this aperiodic 
motion may have been caused by lack of axisymmetry of the particle in question. But it was not possible 
to draw definite conclusions, for several reasons. First the particles were produced by shearing polymer 
microspheres^® which does not produce sufficiently well-controlled shapes. This is a problem because the 
shape of a given particle in the channel cannot be accurately observed. Second it could not be excluded 
that thermal torques affected the orientational dynamics. Third, and most importantly, the setup did not 
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FIG. 1. (a) Schematic picture of the experimental setup, elements are not drawn to scale. Li^ 2 '- lenses in a Keplerian 
telescope configuration. Mi, 2 , 3 : mirrors. Di, 2 : diaphragms. DM: dichroic mirror reflecting the laser beam vertically 
towards the microfluidic system. The microscopic objective focuses the beam on the sample plane. Illumination is 
provided from the top. (b) Coordinate system in the lab frame spanned by orthogonal unit vectors x, y, and z. The 
a;-axis (flow direction) is directed along the channel length, the y-axis along the depth of the channel, and the 2 ;-axis 
along the channel width. The unit vector n points along the major axis of the particle. The polar angle between n 
and the z-axis is denoted by 9. The particles are kept roughly equally far away from the side walls of the channel. 
Since the channel is much wider than deep this means that the y-direction is the shear direction, so that the local 
linearisation of the flow-velocity obeys Ux = sy in the frame co-moving with the centre-of-mass of the particle, where 
s is the shear rate. 


allow to record different orientational trajectories of the same particle. This motivated us to perform the 
experiments summarised in the following, they overcome the problems listed above. 


III. METHODS 

Fig. 1(a) shows a schematic drawing of the experimental setup to observe the orientational motion of small 
particles advected in a microchannel flow. A dilute suspension of microrods in a density-matched fluid is 
introduced into the microchannel. A syringe pump (a standard Harvard Apparatus infuse/withdraw model) 
is used to drive the flow of the particles. The system is placed under an inverted microscope equipped with 
a motorised translation stage. A CCD camera is used to register the orientation of a given particle. The 
initial orientation and position of the particle in the channel is set using an optical trap. 

Fig. 1(b) shows the coordinate system that is used in this paper. The s-axis is parallel to the channel 
length along the flow direction. The y-axis is directed along the depth of the channel. The j/-axis is also the 
optical axis of the microscope objective. The 2 :-axis points along the channel width. The orientation of the 
particle is determined by the unit vector n along the major axis of the particle. 

The experiment is performed with cylindrical glass rods with diameters 3 p.m ± 0.01 |xm (PF-30S, Nippon 
Electric Glass Co., Ltd). The microrods were manufactured as spacers in liquid-crystal devices (PF-30S, 
Nippon Electric Glass Co., Ltd). This requires precise diameters. The lengths of the rods vary between 
approximately 10 qm and 30 qm. An electron-microscope image of the particles is shown in Eig. 2(a). This 
Eigure shows that the end surfaces of the cylindrical rods are randomly inclined and uneven, indicating 
that the particles were obtained by breaking longer glass rods. While this is irrelevant for the intended 
industrial use as spacers, it is important for our application. Random inclinations of the end surfaces break 
axisymmetry: sometimes only very slightly [lower left particle in Fig. 2(a)], sometimes more [c.f. particle 
in the centre of Fig. 2(a)]. We investigate the orientational dynamics of highly asymmetrical particles by 
studying aggregates of glass rods. Following Lewandowski et al. a dilute suspension of microrods is left to 
evaporate in order to produce double particles, Fig. 2(b)-(d). The glass particles have an index of refraction 
of 1.56 and a density of Pp = 2.56 g cm“^. 

To achieve neutral buoyancy the fluid must have the same density as the particles. To this end we mix 
22.2%wt water, 4.4%wt glycerol, and 73.4%wt sodium metatungstate monohydrate (Alfa Aesar GmbH). 
Density matching is achieved by titration while observing the particle under the microscope until the particle 
is buoyant for several minutes. But there are additional sources of error that are difficult to control. For 
example, when the fluid is pumped through the microchannel minute changes in temperature are expected 
to slightly change the fluid density. Such changes could induce small variations in the j/-position of the 
particle that we sometimes observe during the recordings. The suspension is highly diluted in order to avoid 
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FIG. 2. (a) Electron-microscope image of the glass particles used in the experiments. Taken by S. Gustafsson, 

Chalmers, (b-d) Optical-microscope images of the double particles used in the experiments. The angular dynamics 
of the particle in panel (b) is shown in Fig. 5, panel (c) corresponds to Fig. S8, and panel (d) to Fig. S9. The last 
two figures are found in the Supplementary Online Material"*^. 


particle-particle interactions. The mixture has a viscosity of p. = 25mPa s at 20 °C. We estimate the shear 
Reynolds number Reg = pfsa^/p. as follows. The particle length a is of the order of 20 pm. The shear rate s 
is determined by the flow speed Ux which is in turn given by the flow rate 8 pl/min and the cross section of 
the channel, 2.5 mm x 200 pm. This results in m « 0.3mm s“^ and a shear rate of s « 3s“^ at 60pm depth 
assuming a parabolic profile in the j/-direction, and that the channel is much wider than deep. The density 
of the fluid pf equals the particle density, Pp = 2.56 g cm“^. This gives Re^ « 10“^. Inertial effects are thus 
negligible on the time-scale of the experiment. 

The microchannel is produced using standard soft lithography. The process begins by milling a rectangular 
moulding form in aluminum. The surfaces of the mould are mechanically polished. A 10:1 mixture of 
polydimethylsiloxane (PDMS) and sylgard 184 Dow Corning (Sigma Aldrich) is poured on the moulding 
form and allowed to cure for several hours. The PDMS replica obtained after peeling off the mould is sealed 
to a glass slide (thickness 0.17mm) by oxygen plasma bonding. This results in a rectangular channel that is 
40 mm long, 2.5 mm wide, and 150 pm deep. For some measurements, a channel with a depth of 200 pm was 
used, see Table I. The suspension is injected into the channel using thin tubing connected to a syringe pump. 
PTFE tubes (Cole-Parmer) with outer diameters of 0.76 or 1.07 mm are used, the former in connection with 
a plastic syringe (1ml, Terumo), the latter in connection with a glass syringe (500 pi, Hamilton). 

The optical system is built around a Nikon X60 microscope objective (NA 1.0, WD 2mm). The particle 
motion is recorded with a CCD video camera (Leica). The channel is mounted on a translation stage that 
moves the microchannel over the fixed observation microscopic objective. The stage is driven by a stepper 
motor that records the position of the stage. By moving the channel a given particle is kept within the field 
of view of the objective, despite the fact that the particle is advected by the fluid through the microchannel. 

A single-beam optical trap is used to set the initial orientation and position of a given particle. The optical 
trap makes use of the microscope objective (Fig. 1), it provides sufficient magnification to not only visualise 
the particle, but also to trap it with a continuous infrared laser of wavelength 1075 nm (10 W, IPG-Laser 
GmbH). The most efficient way of trapping a glass rod with this setup is to direct the laser beam towards 
one of the two ends of the particle. Different orientations can be imposed on the particle by moving the 
channel sufficiently quickly to cause one end of the particle to leave the trap, yet sufficiently slowly so that 
the other end of the particle is kept trapped. 

The experimental procedure described above ensures stable and reproducible experimental conditions. 
But a number of external factors were difficult to control and could still influence the experimental results 
to a small degree: evaporation due to the presence of unsaturated air in the reservoir, the presence of air 
bubbles in the channel that are difficult to get rid of, and pressure drops ininside the channel due to the 
elasticity of the PDMS structure. 

We use the image-analysis algorithm employed by Einarsson et al. . Images are acquired at a rate of 
100 frames per second. Each frame is stored as 8-bit gray-scale bitmap with 692 x 520 pixels. The pixel size 
is 0.21|.im. For a given frame the image analysis proceeds in three basic steps. First static noise is reduced 
by subtracting the time-averaged intensity from each frame. Then the boundary of the projection of the 
particle into the image plane is detected, and finally an ellipse is fitted to the boundary. Details are given in 
Ref. 27. The output defines the position and the orientation of the particle in the image. The centre-of-mass 
coordinates of the particle in the laboratory frame are determined using the output from the stepper motor 
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recording the position of the stage. 

The main sources of uncertainty in the image analysis are the limited resolution of the camera and 
diffraction. The latter gives rise to a diffuse particle boundary. The uncertainty in the determination of the 
particle orientation is largest when the short end of the rod faces the camera, i.e. for small values of Uz (see 
Fig. 1). 

A typical experiment starts by capturing a particle with the optical trap. The particle is brought into the 
desired location in the x-z-plane by moving the channel. All particles are started close to one of the inlets 
at approximately equal z-distances to both side walls. We verified that the 2 ;-position remains centred, with 
an error typically at most one particle length. This implies that the shear in the z-direction remains very 
small (the channel is much wider than deep). The y-coordinate thus corresponds to the shear direction, 
and the z-coordinate is the vorticity direction. The particle is brought to the desired initial orientation as 
described above and the centre-of-mass of the particle is placed at a depth of 60 mm. The precision in 
determining the initial depth is determine by the depth-of-field of the microscope. It is estimated to be of 
the order of one particle length. Then the particle is released to follow the flow in the microchannel. We 
then invert the pressure gradient so that the particle is advected back in the opposite a;-direction. For each 
orientational trajectory we record both forward and backward dynamics. Since Stokes’ equation is invariant 
under simultaneous pressure inversion and time reversal, the backward dynamics must exactly retrace the 
forward dynamics unless irreversible effects such as inertia or thermal noise affect the dynamics. Examples 
for the resulting video-microscopy recordings of the orientational dynamics in the cc-z-plane can be viewed 
via the MULTIMEDIA VIEW links in Figs. 3 and 5. 

For a given particle this procedure is repeated many times to obtain orientational trajectories with different 
initial orientations. We record the length of the projection of the particle into the a;-z-plane as a function 
of time. We estimate the particle length using the procedure described in Ref. 27. Once the particle length 
is known we can extract the components of the unit vector n determining the orientation of the particle, as 
a function of time. 

We plot the orientation not as a function of time but as a function of distance that the centre-of-mass of 
the particle has traveled through the channel, advected by the flow: 

X{t) = [dt'Uz:{t'). (I) 

Jo 

Here Ux{t) is the instantaneous flow velocity at time t. This transformation simplifies the analysis because 
it accounts for the fact that the shear-rate is time-dependent: the flow velocity changes when the pressure 
is reversed, and in order to avoid inertial effects these reversals must be performed slowly. The invariance of 
Stokes equation under time and pressure reversal implies x(t) = —x{—t). We overlay forward and backwards 
dynamics by plotting the orientation as a function of centre-of-mass position. 


IV. EXPERIMENTAL RESULTS 

Figures 3, 4, and 5 show orientational dynamics of three different particles. For each particle five different 
orientational trajectories are shown, corresponding to different initial orientations [panels (b) to (f)]. In all 
three Figures we show trajectories of the x- and z-components of the unit vector n that points along the 
major axis of the particle. Here Ux is the component in the flow direction, and is the component in the 
vorticity direction. The third component of n is determined by normalisation, jnj = I. 

In panel (a) of each Figure we summarise the orientational dynamics by recording the values of 

U\ 

whenever Ux = 0. We denote the resulting sequence of consecutive riz-values by ni: , i = 1,2,3,.... For 
an axisymmetric particle Jeffery’s equation predicts that = n'z \ shown as the solid line along the 

diagonal in panel (a). 

As mentioned above, the particle is first advected along a stream line of the pressure-driven flow in the 
channel. We then invert the pressure gradient so that the centre-of-mass of the particle is advected back 
to where it came from. For each orientational trajectory we show the ‘forward dynamics’ (blue solid line), 
going from right-to-left in the Figure. After the reversal follows the ‘backward dynamics’ (red dashed line). 
Since Stokes’ equation is invariant under simultaneous pressure inversion and time reversal, the backward 
orientational dynamics must exactly retrace the forward dynamics unless irreversible effects due to inertia 
or thermal noise affect the dynamics noticeably on the time-scale of the experiment. 

Consider first the trajectories shown in Fig. 3, corresponding to particle I. Panels (b) to (f) show orienta¬ 
tional trajectories of Ux and for different initial orientations. In all cases the backward dynamics retraces 
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Position in channel a: [mm] 




FIG. 3. Orientational dynamics of particle 1. (a) Dynamics of riz. Here denote the values of at subsequent 
zero crossings of rix, i = 1, 2, 3,.... The data are taken from panels (b-f). Red o data from panel (b); blue □ data 
from (c); green A data from (d); black + data from (e); magenta * data from (f). Panels (b-f) show orientational 
dynamics as a function of c.o.m.-position x in the channel, Eq. (1). Data in different panels correspond to different 
initial orientations. Solid blue and dashed red lines represent forward and backward trajectories, respectively. The 
flow direction is reversed at a: = 0. The horizontal arrow in panel (b) indicates the period Xp[mm] of the trajectory. 
The video-microscopy recording of the orientational dynamics shown in panel (b) can be viewed at MULTIMEDIA 
VIEW, that of panel (d) at MULTIMEDIA VIEW. 







FIG. 4. Orientational dynamics of particle 2. See caption of Fig. 3 for details. 





























































































































FIG. 5. Orientational dynamics of particle 3. See caption of Fig. 3 for details. The video-microscopy recording of 
the orientational dynamics shown in panel (d) can be viewed via MULTIMEDIA VIEW 


the forward dynamics very well. This shows that neither inertial forces nor rotational diffusion affect the 
orientational dynamics. We attribute the small dephasing visible in each panel to a small density mismatch 
(discussed in Section III), causing the particle to sink (or float), changing the shear rate it experiences. 
Apart from this slight dephasing all orientational trajectories are fairly periodic. For a given na;-trajectory 
the relative variation of the centre-of-mass distance between two consecutive = 0-events (‘half-period’) 
is of the order of 10%. Between different panels we observe variations in the period Xp, between 2.1mm 

and 2.6 mm, caused by the uncertainty in the j/-position mentioned above. Panel (a) indicates that 
is approximately equal to rSz\ 

Fig. 4 shows the orientational dynamics of particle 2 for different initial orientations [panels (b) to (f)]. In 
all cases the backward dynamics retraces the forward dynamics well, at least for a few millimetres. As for 
particle 1 the trajectories show a slight dephasing within each panel. Different panels show quite different 
periods Xp, ranging between 1.7mm and 5.9mm. In panels (d) and (f) we see that the amplitude of 
changes aperiodically. This is also apparent from panel (a) exhibiting a somewhat wider scatter around the 
diagonal than panel (a) in Fig. 3. 

Fig. 5 shows the orientational dynamics of particle 3, an asymmetric double particle. For all initial 
orientations the backward dynamics retraces the forward dynamics fairly well, with the exception of the 
trajectories shown in panel (e). As for particles 1 and 2 we observe a slight dephasing within each panel. 
The dynamics of Ux is quite periodic, but again the periods Xp vary from panel to panel, ranging from 
1.4 mm to 4 mm. In all cases Uz shows distinct aperiodicity. Panel (a) exhibits a larger scatter around the 
diagonal than panel (a) in Fig. 4. 

The orientational dynamics of particles 1, 2, and 3 was obtained under slightly different experimental 
conditions. These are summarised in Table I, as well as the particle properties. Table I also gives information 
about other particles for which we obtained precise orientational dynamics. The corresponding Figures are 
found in the Supplementary Online Material"^^. 

Video-microscopy recordings corresponding to the data shown in Figs. 3(b), (d), and 5(d) can be viewed 
via the MULTIMEDIA VIEW links in the corresponding figure captions. 
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Particle 

number 


Particle 

type“ 

Particle 

length(iim)*' 

Channel 
depth (pm) 

Flow rate 
(pl/min) 

Figure 

Particle 

1 

single 

20.5 

200 

8 

Fig. 3 

Particle 

2 

single 

24.3 

150 

5 

Fig. 4 

Particle 

3 

double^ 

28.7 

200 

8 

Fig. 5 

Particle 

4 

single 

20.1 

150 

7.5 

Fig. 81“^ 

Particle 

5 

single 

29.5 

200 

8 

Fig. 82 

Particle 

6 

single 

17 

150 

5 

Fig. 83 

Particle 

7 

single 

25 

150 

5 

Fig. 84 

Particle 

8 

single 

26.5 

150 

5 

Fig. 85 

Particle 

9 

single 

22.5 

150 

7.5 

Fig. 86 

Particle 

10 

single 

25.5 

150 

5 

Fig. 87 

Particle 

11 

double® 

26 

200 

8 

Fig. 88 

Particle 

12 double-^ 

17 

200 

8 

Fig. 89 


TABLE I. Description of particles and other experimental parameters. “Single or double particle. ^’Particle length 
as extracted from image-analysis algorithm, see Section III. We estimate the error to be of the order of Ipm. “This 
particle is shown in Fig. 2(b). The lengths of the two glass rods are 26.4 and 17.3p.m. ‘^These Figures are in the 
Supplementary Online Material"^^. “This particle is shown in Fig. 2(c). The lengths of the two glass rods are 26 and 
16.5pm. ^This particle is shown in Fig. 2(d). The lengths of the two glass rods are both 16pm. 


V. DISCUSSION 

The results summarised in Figs. 3 to 5 show the orientational dynamics of single and double glass rods. In 
general the particles are not perfectly axisymmetric (as seen in Fig. 2), and therefore do not follow perfect 
Jeffery orbits. In this section we relate our experimental results to the theoretical predictions valid for 
ellipsoidal particles^’^. The particles in our experiment do not satisfy the mirror symmetries assumed in this 
theory, but it is plausible that the effects of breaking axisymmetry predicted by this theory apply at least 
qualitatively to our glass rods. 

The equation of motion for the orientational dynamics of an ellipsoidal particle^can be cast in the form"^^ 


n = On + A (Sn — (n ■ Sn)n) -b 
p = Op + K (E>p — {p ■ Sp)p) + 


---(n • Sp)p, 

KA-l ^ 

( 2 a) 

A(l-iL2) 

KK-l 

( 2 b) 


Following the convention outlined in Section III, n is a unit vector that points along the major axis of the 
ellipsoidal particle. The unit vector p is orthogonal to n, directed along the particle axis corresponding to 
the length b used in the definition of the aspect ratio k (defined in the Introduction). The geometry of the 
ellipsoid is characterised by the two shape parameters A and K that are defined in the Introduction, Section 
I. The matrices S and O are the symmetric and anti-symmetric parts of A, the matrix of fluid-velocity 
gradients. In our case this matrix takes the form 


A = 


0 s 0 
0 0 0 
0 0 0 


( 3 ) 


where s is the shear strength (Fig. I). 

Eqs. (2) are symmetric under the simultaneous exchange of n and p as well as A and K: describing 
the motion of the same particle using a different coordinate system within the particle must result in the 
same dynamics. Note also that the non-linear coupling between n and p involves the strain S only. This 
coupling maintains orthogonality of the two vectors n and p. The anti-symmetric part O just causes a 
solid-body rotation. For axisymmetric particles, AT = 0, so that the tumbling of n becomes independent of 
the dynamics of p (but not vice versa). The resulting equation for n has infinitely many degenerate periodic 
solutions, the Jeffery orbits^. The dynamics of p describes how the particle spins around its symmetry axis. 

When K ^ 0 (and A 7 ^ 0), no general closed-form solutions of Eqs. (2) are known. It is convenient to 
represent the numerical solutions of Eq. (2) in terms of a Poincare surface-of-section^^, recording the locations 
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at which the dynamics intersects a surface in the phase space of Eq. (2). This section is constructed as follows. 
Hinch and Leal ^ have shown that the vector n rotates around the vorticity at a positive angular velocity so 
that one can reduce the dimensionality of the problem by parametrising the orientational dynamics in terms 
of the corresponding angle. A suitable condition^defining the surface-of-section is that n is perpendicular 
to the flow direction, = 0. Following Yarin et al? we take the coordinates in the surface-of-section to be 
'll) and n^, where 'll) is the Euler angle parametrising the spin of the particle around the axis n, and is the 
z-component of n, the cosine of the angle 9 between n and vorticity (Fig. 1). When n^; = 0 we record the 
coordinates To define the Euler angles {6,(l>,'il)) we use the convention of Goldstein^^ and express n 

and p as 



sin 9 sin (p 


— sin xp cos <p — cos 9 sin (p cos xp 

n = 

— sin 9 cos (p 

, and p = 

— sin xp sin (p + cos 9 cos (p cos xp 


cos 9 


cos xp sin 9 


So 6 is the polar angle depicted in Fig. 1. The angle (!> is referred to as the ‘precession angle’. This angle 
measures the direction of the projection of n into the flow-shear plane. Eqs. (2) correspond to Eqs. (2.1) 
and (2.2) in Ref. 2, setting s = — 1 and defining the aspect ratios A and k in terms of the particle axes^ as 
follows: Gx = Xgz and Gy = kGz- 

In the experimental time series shown in Figs. 3 to 5, instances where Ux = 0 correspond to peaks 
in the oscillation in n^. The riz-coordinate in the surface-of-section is therefore easily read off from the 
experimentally observed time series. The angle ip, by contrast, cannot be measured in our experiment 
because we cannot track how the particles spin around n. 

Four different surfaces-of-section are shown in Fig. 6, obtained by numerical integration of Eqs. (2) for 
a large number of different initial orientations, and plotting the sequence of (■)/', riz) evaluated at 

consecutive zero crossings of Ux, labeled by i = 1, 2,3,.... Similar sections are found in Ref. 2. The map 
that gives in terms of is called the Poincare map. 

Fig. 6(a) depicts the orientational dynamics of an axisymmetric particle, K = 0. The coordinate is 
a conserved quantity on the surface-of-section, Jeffery orbits appear as horizontal lines in Fig. 6(a), one- 
parameter families parametrised by xp. In the literature Jeffery orbits are commonly identified by their orbit 
constant C. It is given by the value of tan6* at rix = 0 (see for example Eq. (3) in Ref. I). In this paper 
we characterise Jeffery orbits by Uz = cos0 = 1/(1 -I- G^) on the surface-of-section {'Gx = 0). Fig. 6(a) 
illustrates that the orientational dynamics depends on the initial orientation, determining the value of Uz- 
We remark that the periods of n and xp are not in general commensurate for K = 0. But the tumbling of 
n is independent of that of the angle xp and thus periodic for axisymmetric particles. 

Fig. 6(b) shows the orientational dynamics of a weakly asymmetric ellipsoidal particle [K « 0.095 and 
A = 12/13). We see that Jeffery orbits with riz « ±1 remain almost unchanged. But there are substantial 
changes for smaller values of |nz|, compared with the surface-of-section for K = 0. We see that Uz ceases 
to be a constant of motion, doubly-periodic orientational dynamics results. The most substantial changes 
occur near Uz = 0. The nz = 0-Jeffery orbit is replaced by two fixed points at (0,0) and (±7r/2, 0) on the 
surface-of-section. 

This fact and the surface-of-section patterns in the vicinity of these points follow from the time-reversal 
symmetry of Eqs. (2). The general principle is explained in Section 6.6 of Ref. 44. See also Ref. 46. The 
invariance of Stokes equation referred to in Section III implies that Eqs. (2) are invariant under 


t ^ -t, Ua; -'Gx , and px -Px ■ (5) 

The fixed point (0, 0) is mapped onto itself under this transformation. It follows that the dynamics in 
its immediate neighbourhood can neither be expanding nor contracting. In other words the determinant 
describing the linearised motion in the vicinity of this fixed point. 


det 


dnp+^'> ’ 

dnp'’ 


( 6 ) 


must be of unit modulus despite the fact that the dynamics (2) is dissipative. We find numerically that the 
eigenvalues are exp(±i(j). The point (0,0) is thus an elliptic fixed point, surrounded by a one-parameter 
family of closed orbits that appear as concentric closed curves, much like so-called ‘tori’ in so-called ‘Hamil¬ 
tonian’ systems with area-preserving phase-space dynamics^^. For near-axisymmetric particles Hinch and 
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Leal ^ analysed the corresponding orbits by multiple-scales analysis. As these orbits rotate around the el¬ 
liptic point, the value of Uz changes sign. This doubly-periodic motion may be quasi-periodic or periodic, 
depending on whether the two frequencies are incommensurate or not (corresponding to irrational or rational 
winding numbers, respectively). 

The point (± 7 r/ 2 , 0 ) is a hyperbolic fixed point. We find numerically that the eigenvalues are real and 
opposite, 7 and 7 “^. The fact that the nz = 0 -orbit is destroyed upon infinitesimal perturbation and replaced 
by a discrete set of fixed points of alternating stability is typical for Hamiltonian systems^^. The mechanism 
in our dissipative system is analogous, a consequence of the symmetry (5). 

For larger asymmetries chaotic orientational dynamics occurs, seen as a region with a stochastic scatter of 
points in Fig. 6 (c). Almost all Jeffery orbits are modihed, only those with close to ±1 remain. Whether 
the orientational dynamics is periodic, quasi-periodic, or chaotic depends upon the initial condition on the 
surface-of-section. 

Fig. 6 (d) shows the orientational dynamics for ellipsoidal particles with K = 3/5 and A = 12/13. This 
value of K is similar the corresponding value for the double particle 3 (Figs. 2(b) and 5), but this particle 
has A « 40/41. We do not seek quantitative correspondence between these parameters because the particles 
used in our experiment are not ellipsoidal. This prevents us from drawing quantitative conclusions, but still 
allows for a qualitative comparison between theory and experiment. The orientational dynamics displayed 
in Fig. 6 (d) is either on tori or chaotic. Motion on tori can occur with large amplitudes, so that changes 
from Uz ~ —1 to Uz ~ 1. 

How are these observations reflected in the experimental time series shown in Figs. 3 to 5? Since we 
cannot measure the angle xp in our experiments, we concentrate on the dynamics of n^. Fig. 7 shows the 
range of changes of Uz in one iteration of the Poincare map for particles with different degrees of asymmetry, 
determined by numerically recording the changes n/ —nk along orientational trajectories. The diagonal 
corresponds to symmetric particles where The larger the asymmetry, the larger is the range 

of — n^z^ that may occur, reflecting the properties of the surfaces-of-section shown in Fig. 6 . 

Fig. 7 can be directly compared with Figs. 3(a) to 5(a). Our experimental results show that the Poincare 
map may scatter significantly around the diagonal. The above discussion explains that this is a consequence 
of doubly-periodic and possibly chaotic orientational dynamics of asymmetric particles. The range of scatter 
differs between different particles, a consequence of different degrees of asymmetry. We see that the scatter 
is largest for the double particle 3, with K « 3/5 [Fig. 5(a)]. Particle 1, by contrast, shows only negligible 
scatter around the diagonal [Fig. 3(a)]. We infer that this particle is highly symmetric, K is very small. 
The data shown in Fig. 3 are consistent with the conclusion that particle 1 follows Jeffery orbits, rZz remains 
approximately constant in the surface-of-section. But we remark that the shape of the nz ~ 0-orbit does not 
look like a Jeffery orbit for an axisymmetric particle. We cannot exclude that this is due to possible non- 
ellipsoidal deviations from axisymmetry. The surface-of-section dynamics is most sensitive to such shape 
perturbations near = 0. A more likely explanation is that the shape of the trajectory is a consequence 
of systematic (and reproducible) tracking errors due to diffraction and finite pixel size. These errors are 
largest for small values of n^, as discussed in Section HI. 

Particle 2 is also a single glass rod, but it shows a somewhat larger range of scatter around the diagonal 
[Fig. 4(a)]. We attribute this to a more substantial breaking of axisymmetry at the tips of the particle [as seen 
for instance in Fig. 2(a)]. Particle 2 shows fairly periodic motion for riz « 1 and distinctly doubly-periodic 
orientational dynamics for small values of nz- This confirms that the initial orientation determines whether 
the particle tumbles periodically or aperiodically. Note also that Uz changes sign along the trajectories that 
remain near = 0 [Fig. 4(b), (d), and (f)]. All of these observations are qualitatively consistent with the 
surfaces-of-section shown in Fig. 6 (b) and 6 (c). 

For particle 3 the value of changes sign for all orientational trajectories shown. This is consistent 
with the fact that the Fig. 6 (d) shows predominantly this type of motion. Fig. 5(e) is consistent with 
chaotic orientational motion in the stochastic layer around the elliptic island. The surfaces-of-section in 
Fig. 6 show that there are two types of doubly-periodic motion: either nz has the same sign, or its sign 
changes periodically. The trajectory in Fig. 5(e) exhibits both behaviours, indicating chaotic dynamics. 
Now consider the trajectory shown in Fig. 5(c). It is not periodic or doubly-periodic, but the reversal works 
well (at least initially). But we cannot conclude that the dynamics is consistent with surfaces-of-section 
discussed above. The riz-values at = 0 in Fig. 5(c) change from approximately 0.25 to 0.75 in modulus. 
Explaining this behaviour in terms of chaotic dynamics on the surfaces-of-section requires a large stochastic 
region, larger than the one shown in Fig. 6 (d). This might mean that the trajectory shown in panel (d) 
is a doubly-periodic piece of a chaotic trajectory that may show different behaviours at larger times. But 
to determine whether these behaviours can be explained by chaotic dynamics would require to derive and 



FIG. 6. Poincare surfaces-of-section for A = 12/13 and different values oi k [K = {k? — 1)/(k^ + 1)]. The angle 
is defined np to {tp + 7r/2)niod7r. 



FIG. 7. Range of Poincare map for A = 12/13 and different values of k corresponding to the values used in Fig. 6, 
K = + 1). The data shown was obtained by by iterating the Poincare map once for randomly chosen 

initial conditions. The initial conditions were determined by sampling ip uniformly over the surface-of-section for 
given values of Uz- 
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numerically integrate the orientational equations of motion for the precise shape of the particle, and to 
experimentally record the angle •0. Since the current experimental setup does not allow to reliably extract 
how the angle if) changes, we plan to perform experiments with small triangular platelets that will allow to 
record the angle •0. 

The winding number of the trajectory shown in Fig. 5(d) is roughly 7 corresponding to seven Ux = 0- 
crossings while the trajectory winds approximately once round the elliptic fixed point in the centre of the 
surface-of-section. The winding number of the trajectory in Fig. 6(e) is roughly 4 corresponding to four 
Ux = 0-crossings while the surface winds once around the elliptic fixed point. These observations are in 
qualitative agreement with the fact that the winding numbers of tori winding around the elliptic fixed points 
increase as the distance from that point increases (corresponding to larger maximal values of Uz). 

Additional results for seven more particles are shown in the Supplementary Online Material"^^, Figs. SI 
to S7. In general the results shown in these Figures support the observations and qualitative conclusions 
summarised above. But the trajectory shown in Fig. S7 is difficult to reconcile with the surfaces-of-section 
shown in Fig. 6. Fig. S7(c) shows an orbit near Uz = 0 where Uz appears not to change sign. But panels (d) 
and (f) in Fig. S7 show distinct sign changes, not consistent with the surfaces-of-section shown in Fig. 6. 

It was pointed out in Section IV that the periods Xp observed in the orientational trajectories can differ 
substantially between different trajectories of the same particle. Consider for instance particle 1 (Fig. 3). 
For an axisymmetric ellipsoidal particle {K = 0) the Jeffery period (time units) is given by^ 


Tp = TT 


-f 1 
sX 


( 7 ) 


For K = 0 this is twice the return time to the surface-of-section. For K 0 the return time depends upon 
the starting position on the surface-of-section, but for small values of K the deviations from (7) are small, 
of the order of K. For the nearly axisymmetric particle 1 we can use Eq. (7) to estimate the period Xp in 
Fig. 3. Using the parameters from Table I and assuming that the particle was located at a depth of 60 pm, 
Eq. (7) gives Xp « 2.3 mm. This is consistent with the range of periods observed in Fig. 3(b) to (f), namely 
2.1 mm to 2.6 mm. We infer that the precision in determining the depth at which particle 1 moves through 
the channel is of the order of one particle length. The variations in periods observed in Figs. 4 and 5 indicate 
that the actual depths vary more between different panels than in Fig. 3. 


VI. CONCLUSIONS 

Theory and numerical simulations predict that the orientational dynamics of small neutrally buoyant 
particles in a shear flow is very sensitive to breaking of axisymmetry^’^’^^. Axisymmetric particles tumble 
periodically on Jeffery orbits, but when the symmetry is broken, the Jeffery orbits are modified. Depending 
on the initial orientation, periodic, doubly-periodic, or chaotic tumbling may result. 

In order to experimentally verify these predictions it is necessary to use particles of well-defined shape, 
to ensure that inertial effects and rotational diffusion are negligible, and to compare different orientational 
trajectories (corresponding to different initial orientations) of the same particle. 

In this paper we describe measurements of the orientational motion of small glass rods suspended in a 
micro-channel shear flow. The measurements are precise enough to allow for a definite comparison with 
theory^’^. First, the glass rods have highly symmetric circular cross sections. Slight imperfections at the 
ends break the axisymmetry weakly for some particles [Fig. 2(a)]. We also use strongly triaxial particles 
designed for our purpose by joining two glass rods. Second, by reverting the pressure-driven flow we observe 
that the orientation retraces its trajectory for many periods. This means that neither inertial effects nor 
rotational diffusion matter on this time scale. Third by means of an optical trap we change the orientation 
of the particle, making it possible to measure different orientational trajectories for the same particle. 

Our results confirm the theoretical predictions^’^. We observe periodic and doubly-periodic tumbling, and 
how the nature of the orientational dynamics depends upon the initial orientation, and on the degree to 
which axisymmetry is broken. Our measurements are consistent with the features of the surfaces-of-section 
shown in Figs. 6 and 7. 

In the future we plan to experimentally map out Poincare surfaces-of-section (Fig. 6). With the present 
particles this is not possible because we cannot resolve how the particle spins around its major axis. We 
therefore plan to perform corresponding experiments using micron-sized triangular platelets. We expect 
that it will be possible to resolve the tumbling and spinning dynamics of these particles using the methods 
described in this paper. 
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FIG. SI. Orientational dynamics of particle 4. (a) Dynamics of riz. Here nz(i) denote the values of riz at 

subsequent zero crossings of rix, i = 1,2,3,.... The data are taken from panels (b-f). Red o data from panel (b); 
blue □ data from (c); green A data from (d); black + data from (e); magenta * data from (f). Panels (b-f) show 
orientational dynamics as a function of c.o.m.-position x in the channel, Eq. (1) in the main text. Data in different 
panels correspond to different initial orientations. Solid blue and dashed red lines represent forward and backward 
trajectories, respectively. The flow direction is reversed at a; = 0. 




FIG. S2. Orientational dynamics of particle 5. 
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FIG. S3. Orientational dynamics of particle 6. 
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FIG. S4. Orientational dynamics of particle 7. 
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FIG. S5. Orientational dynamics of particle 8. 
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FIG. S6. Orientational dynamics of particle 9. 










































































































































5 





Position in channel ®[mm] 



TT^ 

(e) 


n 

1 

1 

1 

n 

1 

1 


n 

1 


1 


1 

L 

/ , 

L 

1 

1 

ii 

1 

1 

1 

ii 

y=i 

1 

[j 




FIG. S7. Orientational dynamics of particle 10. 
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FIG. S8. Orientational dynamics of particle 11. 
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Position in channel a: [mm] 

FIG. S9. Orientational dynamics of particle 12. 















